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Abstract 

We reconsider the Feynman-Kishnger-Ravndal model apphed to neutrino-exci- 
tation of baryon resonances. The effects of lepton mass are included, using the formal- 
ism of Kuzmin, Lyubushkin and Naumov. In addition we take account of the pion-pole 
contribution to the hadronic axial vector current. Application of this new formalism 
to the reaction i^f^+p ^ + A"*""*" a.t Ej, ^ 1 GeV gives a suppressed cross section at 
small angles, in agreement with the screening correction in Adler's forward scattering 
theorem. Application to the process + P ^ t~ + A+"*" at E'y ~ 7 GeV leads to the 
prediction of right-handed t~ polarization for forward-going leptons, in line with a 
calculation based on an isobar model. Our formalism represents an improved version 
of the Rcin-Sehgal model, incorporating lepton mass effects in a manner consistent 
with PCAC. 



1 Introduction 

A new generation of neutrino experiments is under way that is exploring low energy neu- 
trino reactions such as i^^n — > n~p and v^p — > /i~p7r+ with unprecedented statistics and 
detectors of high resolution [1]. As an example, the MiniBoone experiment has reported 
preliminary results on 40000 events of the type f^p fi~pTr~^ [5]. These data are providing 
incisive tests of theoretical models and revealing features that pose new challenges to the- 
ory. One such feature has been referred to as a "deficit of forward muons" in the reaction 
I'fi+ip, n) fi~ + {p, n)-|-7r+ and its coherent counterpart z/^-|-Nucleus /i~-|-7r+-|-Nucleus 
[3ll3|. In a recent paper [5] it was pointed out that the inclusion of a non-zero muon mass 
causes a suppression of the coherent process in the forward direction, as a consequence of 
a destructive interference induced by spin-zero pion exchange. 

On a different front, lepton mass effects are of obvious importance in reactions like 
VrP T~p-K^ which are being discussed in connection with experiments to detect — > i^r 
oscillations of atmospheric or laboratory neutrinos [6] . The large mass of the r lepton has 
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implications for the angle and momentum distribution of the r , as well as its polariza- 
tion. Thus a formalism that includes non-zero lepton mass is indispensable. 

In this paper, we propose an extension of the Rein-Sehgal (RS) model [7] for resonance 
production vi + N ^ /" + A^* that takes account of lepton mass effects. This model, based 
on the original work of Feynman, Kislinger and Ravndal [H [9], has been successful in 
describing data in a variety of neutrino experiments over the past 25 years, and has served 
as a code for simulating single pion production in the resonance region up to ~ 2 GeV. 
The attractive feature of the model is its economy: using as input the vector and axial 
vector form factors of the quasi-elastic channel v^n — > /U~p, it provides a unified descrip- 
tion of resonance production, embracing nearly 20 resonances in the above mass range. 
The predictive power of the model derives entirely from SU(6) spin-flavour symmetry and 
the conservation properties of hadronic weak currents. 

In introducing lepton mass corrections, we will make use of a formalism developed 
by Kuzmin, Lyubushkin and Naumov [TU], but will modify it in a significant way that 
respects the PCAC property of the weak axial vector current. 



2 Brief Review of the RS Formalism [mi = 0) 



We recall the basic features of the RS model in the limit mi = [TJ [U [9]. The matrix 
element for a typical resonance excitation process + iV — > ;U~ -|- N* is written as 



T{v^N ^ ^l-N*) = ^ [W7^(l - 75)«.j (iV*| J+(0)|iV) 



where G = Gp cos 6^, and the hadronic current operator is expressed as 
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M denoting the resonance mass. The lepton current is expanded in the rest frame of the 
resonance (RRF) as 



RRF 
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f -|- V 2uve{ 
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where = —q^ is the momentum transfer, and u and v are kinematical factors depending 
on the initial and final lepton momenta (see Eq.(2.9) of [^). The vectors e^, e^, may 
be regarded as polarization vectors of the virtual intermediate boson, corresponding to 
left-handed, right-handed and scalar polarization, 



e^ = -=(0,1,-^,0), e^ = -=(0,-1,-^,0), e 
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where = (z^*,0,0,Q*) is the momentum transfer 4- vector in the resonance rest frame. 
The matrix element then takes the form 



T{„^N ^ fi"N*) = -AGUE, I J^{N*\uF. - vF+\N) + ^V2^{N*\Fo\N) i (5) 



with 
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24F, = Ft + ^F, 



Q* 

(In Eq.Q, Mn denotes the nucleon mass, and q the 3-momentum transfer in the lab.) 
The resulting differential cross section is 



with 



3z 



(8) 



K denoting the conventional "flux factor" k = iW"^ — M-^)/2Mj\f. The helicity amplitudes 
/±[2,,| = {N,j, ± l\F±\N*,j,) and /o± = {N,±^\Fo\N* ,±^) are listed in Table II of [7] 
for all known resonances up to = 2 GeV. They are expressible in terms of 3 functions 
proportional to {Q"^) (these are called T^, and S) and 4 functions proportional to 
G^{Q'^) (these are called T^, R^, B and C). These 7 functions (called reduced matrix 
elements or dynamical form factors) are listed in Eq.(3.11) of Ref.[7]. 



3 Lepton Mass Effects due to Nonconservation of Lepton 
Current 

The RS formalism valid for massless muons has been extended by Kuzmin, Lyubushkin 
and Naumov (KLN) to the case mi ^ [10]. The first difference is that the lepton in the 
final state can now have helicity + or — , since the lepton current is no longer conserved. So 
there are six helicity cross sections (^^^\ ^^s^i A = + or — . In addition, the kinematical 
coefficients multiplying these cross sections are no longer identical to the factors u^, v^, 
2uv that occur in the massless case (Eq.d?])), since the massive lepton satisfies a different 
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energy-momentum relation. In the notation of KLN, the components of the lepton current 
in the RRF depend on A, and may be written as 
1 



^V) =Aa)1^V1-Acos0(M^ -El- XPi) 



Jx{X) 
Jz{X) 

where 



W 
1 



^(A)|^Vl + ACOS0(P; - XE^) 

iXA(^x) Vl + A cos 9 

Aa)^777V1- Acos0 [{E^ + XPi){Mn - Ei) + Pi{XE^ + 2E^ cos 9- Pi)] 



(9) 



\q\W 



A, 



!(;,) = ^E,{Ei - XPi), (10) 

and the symbols E,^, Ei, Pi and 9 are variables in the lab frame. This modified lepton 
current can be expanded, as before, in terms of three polarization vectors corresponding 
to left-handed, right-handed and scalar polarization. 
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where 



eg =-^(0,1,-^,0) 
4 =4f (0,-1,-^,0) 



(12) 



Q(A),o,o,z., 
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Note that the vectors eg and e^ are the same as in the case mi = 0. The novelty is 
in the scalar polarization e"^^ whose components depend on the lepton helicity A. The 



coefficients c'^\ c'd^ and cL'^^ are given by 
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while the components Q*(^x) ^(A) ^'■^ 



(A) 



(Note that starred quantities refer to components in the RRF.) With these definitions, the 
differential cross section is 
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0.) Remarkably, the helicity cross sections f)^^^ can be calculated exactly as 
in the RS model (using Eq.(2.15) and Table II of Ref.[7j) provided three of the dynamical 
form factors are modified as follows: 



S Skln ■ 
Bkln ■ 



^(A)^ 



Q 



(A) I 



1 + 



Q2 3W\G^{Q 



Mn J 6|q| 



B 
C 



? ( Q(A) + ^1 



Ckln ■ 



1 U* 

3 ^ aMN 



(16) 



(A) 



aMN ' 3aMN ) 



+ 



Here, 



E*. 



W 



1 + 



ZG^ (Q^) 
2W\q*\ 



TV 



(17) 



2MnW 

As in Ref.[7], Z ~ 3/4 is the renormalization factor for the axial vector current, = 1.05 
GeV^ is the slope of the baryon trajectory and n is the number of oscillator quanta in 
the final resonance. The form factors G^ (Q^) ^"^^ {Q"^) taken to be dipoles 



G^ (Q2) = ^1 + 
GeV, Ma = 0.95 GeV. 



G^ (Q2) 



1 + 



, the default values being My = 0.84 



4 Lepton Mass Effects due to Pion-Pole Term in Hadronic 
Axial Vector Current 

The KLN formalism described in SecjS] takes account of lepton mass corrections stemming 
from the modification in the components of the lepton current, keeping the hadronic 
current unchanged. These corrections can be implemented within the RS formalism by 
redefining the dynamical form factors S, B and C, as indicated in Eq. (ll6p . There is, 
however, a further effect associated with the nonvanishing lepton mass that remains to be 
considered. As discussed by Ravndal [9j, the axial hadronic current has, besides the quark 
current defined in Eq.(4.6) of that reference, also a pion-pole contribution, dictated by 
PCAC, which modifies the axial current as follows: 

A,^A; = A, + q, J^^^ (18) 
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As long as the lepton mass is neglected, the additional term in ^4^ multiplied by the 
lepton current gives zero, and so has no effect. If ^ 0, however, the pion-pole term 
does contribute. In particular the divergence of is 



As shown by Ravndal [9], this modification of the axial vector current leads to a matrix 
element for the quasi-elastic process + n ^ /i~ + p of the form 



with F2^(Q2) = ^A(g2 



+ (20) 
ml + Q"^ 



Neglecting the small correction m^/M]\], this result implies an induced pseudoscalar form 
factor 

i^2^(0) = ^5a(0), 5a(0) = 1.25 (21) 



which agrees with the PCAC result [-^2^(0)] p^^^^^ = V^gNNTrf-n/m^ provided (/^(0)2mAr = 
V^gNNirf-K, which is the Goldberger-Treiman relation {q^j^t^/^t^ ~ 14, /tt = 130 MeV). 
Thus the inclusion of the pion-pole term in the axial current is mandatory for a satisfac- 
tory description of the quasi-elastic matrix element, and will therefore affect the amplitudes 
for resonance-excitation as well. 

The effects of the pion-pole term can be incorporated by recalculating e^A^^ and q^A^^ 
in the manner described by Ravndal (see Eqs.(4.7) and (4.8) of [9]), using the modified 
axial vector current, Eq.p^. Only the reduced amplitudes B{Q'^) and C{Q^) are affected, 
and the corrected form, relative to that computed by Kuzmin et al., is 



r(a) _„(a) ZGa (Q') / ^}hv_^N 



Mat a 



BRS-^KLN+ 2WQ* y^^^^" -"W^ ) m2+Q2 



(22) 



Thus the complete set of dynamical form factors defining the extension of the RS model 
to take account of muon mass effects is T^, , T"^, R"^ (these are unchanged) together 
with S^^Ij^, B^^^g and C^^^g given in Eqs. (fT6]) and ([22|) . 

5 Application to z/^ + j9 ^ /i + 

At energies of order 1 GeV, the channel /i^pvr"'" is dominated by the 

A''~''~(1232) resonance, and serves as an interesting testing ground for the model discussed 
in this paper. 
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Figure 1: Differential cross section da/dco&O^ in units of 10"'^^ cm^ for the reaction 
v^p ii~TT^p aX E = 0.7 GeV. a) scalar cross section for the three models 
discussed in the text, b) transverse cross section. Here BRS and KLN give 
identical results. 



5.1 Angular Distribution da/dcos9^: Transverse and Scalar Contribu- 
tions 

To discuss the angular distribution of the muon it is useful to decompose the cross section 
da/dQ'^dW^ into transverse and scalar pieces, 



where 




Examination of the kinematical coefficients , Cj^ and given in Eq. ()13p shows that 
for forward scattering {6^ = 0), the only coefficients that survive, are and \ The 
coefficient represents a helicity-flip amplitude, which means, that the only surviving 
contribution to {da /d cos 0^)^^ at 0^ = is a helicity-flip cross section which vanishes when 
0. By contrast, the scalar contribution (dcr/d cos 0^)^^ at 0^ = is determined by 

Cg \ and the corresponding helicity-conserving cross section o"^ \ which involves the dy- 
namical form factor C^^RS (^^) given in Eq. ()22p . This form factor is particularly sensitive 
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to the lepton mass because of the pion-pole contribution. 



The scalar contribution to da/dcos9^ for Ei, = 0.7 GeV is plotted in Fig. la. The 
curve labelled 'RS' is based on the model of Ref.[7], with matrix elements calculated for 
massless muons, but muon mass included in the phase space. The curve 'KLN' describes 
the result of Ref.fTHj, based on the dynamical form factor Ck^ln{Q'^)^ which contains 
muon mass effects excluding the pion-pole. The full correction is contained in the curve 
labelled 'BRS' based on the dynamical form factor Cbrs {Q'^) calculated in the present 
paper. Notice that the scalar cross section at small angles is suppressed by the effects of 

7^ (compare 'BRS' with 'RS'). 
The transverse cross section {da/dcos9^)^^ is plotted in Fig.lb, where we compare the 
'RS' model with the lepton mass corrected result (here there is no distinction between 
'KLN' and 'BRS'). The difference relative to 'RS' reflects the non-conservation of lepton 
helicity, the cross section containing both A = + and A = — . The non- vanishing result in 
the forward direction represents the helicity-flip component c"i\o = 0) 7^ 0. 

Finally, in Fig. 2, we add the scalar and transverse pieces to show the full angular 
distribution da /d cos 9^. The muon mass effect is contained in the difference between 
'BRS' and 'RS', and amounts to a suppression in the angular interval 0.9 < cos^^ < 1.0 
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Figure 3: Differential cross section per nucleon da/dcosO^ in units of 10^'^® cm^ for the 
reaction + ^^C + tt^ + ^^C. The sohd Une represents the coherent 

Rein-Sehgal model as calculated from ^7} . The dashed line includes the Adler 
screening factor a) = 0.7 GeV b) = 1.3 GeV. 



of 16.5%. We have verified that in the limit in which the only mass-corrections retained 
are those in which a factor is accompanied by a pole-term ^i^qi , the cross section in 
the forward direction is 



'da{e = oy 


BRS ^ 


d cos 6^ 





Q 



2 

min 



2QI 



+ mi 



^da{9 



0) 



dcos 6n 
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where Q 



2 

min 



m'iy/{l — y), y = v/E^, and v"^ » Qmin- have thus recovered the 



"screening" factor (l — ^Qmin/iQmin + '^^)) which appears in Adler's forward-scattering 
theorem |llj , and which was invoked in Ref . [5] in order to explain the forward muon deficit 
observed in Ref. [3] . 



5.2 Comparison of Coherent and Incoherent Reactions 

The angular distribution of the n~ produced in the incoherent process i'^ + p ^ fJ,~ + 
A"*""*" — > /i~ -|- vr"'" -|- p is markedly different from that of muons produced in the coherent 
process 

z/^ -I- Nucleus fi~ +TT~^ + Nucleus (26) 

The dynamics of the coherent reaction is governed by the divergence of the weak axial 
vector current, and the muon angular distribution may be obtained by appealing to Adler's 
PCAC theorem for small angle scattering. An explicit model for this process was developed 
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in Ref.|12j and leads to the differential cross section (for = 0) 



dQ^dydt 2tt^-"' y Wtt 



{l + r^)( rln'^ ^ e~'^\'\FaUE. = E,y) (27) 



where and y are the momentum and energy transfer variables at the lepton vertex, 
t is the square of the momentum transfer at the nucleus, and the various factors are 
defined in [12]. Using the above formula we have calculated da /d cos 9i for the reaction 
+ ^^C ^ /i~ + 7r+ + ^^C at El, = 0.7 GeV, using empirical data on the total and elastic 



I', 



TT^p cross section to evaluate atot{'^'^ J^) and Fabs, with M denoting an average nucleon. 
The resulting angular distribution da /d cos 9 fj_ is shown in Fig. 3a. When the effects of 
muon mass are included, the differential cross section in Eq. ()27p is multiplied by the Adler 
screening factor [H [H] 

This leads to a damping effect in da /d cos 9^ at small angles, exhibited in Fig. 3a. The 
corresponding effect at E^, = 1.3 GeV is shown in Fig. 3b. The contrast between the highly 
peaked distribution in Fig. 3a and the corresponding broad distribution of the incoherent 
process plotted in Fig. 2 suggests that the variable cos0^ is an appropriate choice for 
distinguishing between the two processes. 



5.3 Distribution in Q'^ 

The variable has also been used as a probe of the form factors involved in the dynamics 
of A"^"*" production, and the distribution in this variable is shown in Figs. 4a, 4b for the 
energies Ei^ = 0.7 and 1.3 GeV, using the BRS model. Also shown is the combined (co- 
herent + incoherent) distribution, where the coherent cross section includes the screening 
correction, Eq. (j28p . An estimate of the suppression due to muon mass effects may be 
obtained from Table 1. In the interval < 0.1 GeV^, the incoherent cross section for 
El, = 0.7 GeV is suppressed by 14%. For the coherent process, nearly 80% of the cross 
section is in this Q^-bin, and the suppression due to screening is 16%. Note that this 
suppression refers to a comparison of the BRS calculation with an unscreened RS model, 
in which a non-zero muon mass is retained in the phase space. A stronger suppression is 
obtained in Ref.[5], where the screened result is compared with an RS calculation with 
= everywhere. 

It may be added that for neutrino reactions in nuclear targets, a further source of 
suppression at low is a possible "Pauli-blocking" effect, discussed, for example, in the 
papers in Ref.|13j. Estimates based on a Fermi gas model indicate that the incoherent pro- 
cess could undergo an additional suppression of ^ 5% at < 0.1 GeV^, as a consequence 
of such nuclear effects. 
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Figure 4: Differential cross section da/dQ^ in units of 10^'^'^cm^/GeV^ for 7r+ produc- 
tion by muon neutrinos in charged current reactions, a) i?^ = 0.7 GcV b) 
En = 1.3 GeV. Tlie incoherent cross section for an average nucleon was cal- 
culated in the A-approximation, using (J^g/^ — (5/9)(tJ^q/. 
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Coherent scattering 
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E = 1.3 GeV 
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CrQ2<0.1 


a 
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C^Q2<0.1 
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0"Q2<o.l 


RS 


0.227 


0.049 


0.504 


0.073 


RSC 


0.173 


0.141 


0.305 


0.242 


BRS 


0.194 


0.042 


0.483 


0.067 


RSA 


0.147 


0.118 


0.283 


0.222 



Table 1: Integrated cross section for incoherent and coherent neutrino reactions in units 
of 10~'^*cni^. The index < 0.1 indicates integration of dajdQ^ between 
= and = 0.1 GeV^. The coherent cross section per nucleus has been 
calculated for ^^C, taking a Carbon radius of 2.42 fm. The line labelled 'RSC' 
refers to the coherent Rein-Sehgal model as calculated from ^7\i . The line 
labelled 'RSA' includes the Adler screening factor (j28p . 
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6 Application to + p ^ r + A 




1 2 3 4 5 6 7 



Figure 5: The absolute value of the t" polarization at 0° versus r momentum pr for 
three different models used in calulating the cross section of VtP ^ r^pn'^. 
The solid line (labelled BRS) refers to the model presented in this paper, 
whereas the dashed line (labelled KLN) refers to [TU] and the dashed-dotted 
line (labelled HMY) to [E] • The neutrino energy is = 7 GeV. 

The consequences of lepton mass corrections are quite dramatic for resonance produc- 
tion by I'r- In particular, it was noted by Kuzmin et al. [H] that in the KLN version 
of the model, there are significant dynamical effects of ^ 0, over and above the 
purely kinematical (phase space) effects. Distributions studied in [Hj include da/dQ'^ 
for ly-rP — > T~pTT~^ for various neutrino energies, and da /dp-j-d cos Oj- at various neutrino 
energies and angles Or- One can expect that the 'BRS' version of the model will yield 
predictions that differ from those obtained in [T3]. In one instance, Ref.[14j found a result 
that was in disagreement with a calculation based on an isobar model [15] for the process 
i/^ -\- p ^ T~ + A"'"^, which uses a Rarita-Schwinger formalism, and phenomenological 
form factors constrained by CVC and PCAC. The observable studied was the degree of 
polarization of the r~ in z^,- + p ^ r" + A"^"*", defined as the fractional difference between 
helicity A = + and A = — of the final t~ lepton. The configuration chosen was 6r = 
(forward production) and the polarization calculated as a function oipr, the r momentum, 
at an energy = 7 GeV. 

We have investigated this discrepancy using the present (BRS) version of the RS model 
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which includes lepton mass corrections due to the pion-pole. One expects that the dynam- 
ical form factor Cbrs would have an impact for forward t~ production (9^ = 0), where 
the screening effect of the pion-pole term should be significant, and the reaction could be 
dominated by the helicity-flip term o"j-^\ In Fig. 5, we show our results for the polarization 
of T~ , defined as 

dajX = +)- dajX = -) 
^ ~ da{X = +) + da{\ = -) 

as a function of the r momentum p^, for 0^ = and Ey = 1 GeV. The curves compare 
the polarization obtained in the isobar model with the KLN and BRS versions of the 
relativistic quark model. The BRS result for is positive (right-handed helicity), and 
compatible with the isobar model, whereas the KLN result changes sign. (The cusp in 
the KLN curve is due to the fact that Fig. 5 plots the modulus of the polarization.) We 
see this as an indication that the incorporation of the pion-pole into the RS model, in 
accordance with the prescription of Ravndal, produces lepton mass effects which are in 
accord with phenomenological approaches based on the conservation properties of weak 
vector and axial vector currents. 
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